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1 . INTRODUCTION 

The positional accuracy of photo-triangulated ground points can 
be evaluated either by the use of accurate check points or by statisti- 
cal analysis based on the principles of error propagation. In conven- 
tional engineering applications of aerotriangulation, some known 
ground points are usually withheld from the aerotriangulation solution 
and are used to serve as ultimate check on the accuracy of the aero- 
triangulation process. Complete reliance on statistical analysis 
procedure is not an advisable or corranon practice. 

However, in lunar phototriangulation, there is a complete lack 
of accurate ground control points. The accuracy analysis of the 
results of lunar phototri angulation must, therefore, be completely 
dependent on statistical procedure. It was the objective of this 
investigation to examine the validity of the commonly used statistical 
procedures, and to develop both mathematical techniques and computer 
softwares for evaluating 1) the accuracy of lunar phototriangulation; 

2) the contribution of the different types of photo support data on 
the accuracy of lunar phototriangulation; 3) accuracy of absolute 
orientation as a function of the accuracy and distribution of both 
the ground and model points; and 4) the relative slope accuracy be- 
tween any triangulated pass points. 

These statistical techniques were applied to evaluating the tri- 
angulation and mapping accuracy of the mapping system flown in the 
Apollo 15, 16 and 17 missions. The system included a 76-mm. metric 
camera with a 74®x74° field of view and a frame size of 114x114 mm. 

The camera had a glass reseau plate containing 121 reseau crosses at 
the focal plane. The camera provided 1:1,450,000 photography at a 
nominal altitude of 110 km. with a ground resolution of 20 m. A 
76-mm. stellar camera took pictures of the star field to provide 
attitude controls for the mapping camera. The two cameras were timed 
for simultaneous exposures. Scale control was provided by a laser 
altimeter which had a designed range accuracy of +2 m. Also included 
in the system was a 610-mm. optical bar panoramic camera. It had a 
10'’46'xl08‘’ field with an image size of 114x1,148 mm. At an altitude of 
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no km, its photography had a scale of approximately 1:180,000 at the 
center of the frame with a resolution of about 1 to 2 m. The panoramic 
camera was tilted alternately forward and backward 12.5° so that 
consecutive frames of similar tilt had 10% overlap and stereopairs 
had 100% overlap. It provided sterocoverage of a strip approximately 
330 km, wide, centered on the ground tracks. (5) 

The 76-mm. mapping camera photography is being used by NASA to 
establish a network of lunar controls by means of phototriangulation. 
The high-resolution panoramic photography is used to produce large 
scale lunar maps. The pass points obtained from the triangulation of 
the 76-mm. photography will be used as ground controls in the absolute 
orientation of the stereomodels from the panoramic photography. In 
this investigation, the relative accuracy of the triangulated pass 
points was evaluated by a simulation process; and the accuracy of the 
absolute orientation process was evaluated using real data provided 
by NASA. 

Mr. Nick 6. Yacoumelos was the research assistant in this project 
when he was studying for the doctorate degree at the University of 
Illinois at Urbana-Champaign. He did an extensive literature research 
on the method of propagation of variance, and assisted ably in both 
data processing and data analysis. He is presently an Assistant 
Professor at the Lowell Institute of Technology. 

Mr. Robert Hill of the Mapping Sciences Branch at the Lyndon B. 
Johnson Space Center was the technical contract monitor for NASA, 
AS-ll-Bl data on one model of 76-mm. photography and three models of 
panoramic photography from the Apollo 15 mission was provided by the 
Lockheed Electronics Company in cooperation with Mr. Robert Hill. 
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2. THEORETICAL BASIS OF THE METHOD OF PROPAGATION OF VARIANCE 
In the application of the least-squares method in analytical 
photogrammetry, the variance-covariance matrix of the unknown param- 
eters is commonly computed as the inverse of the normal equation 
multiplied by the standard error of unit weight. For example, in 
the simultaneous solution for phototriangulation, the solution 
equations may be sunmarized as follows (4): 

6 = N"^C; (2.1) 

and 

6j = Nj’^Cj - for all pass points j; (2.2) 

where 6 is a matrix of the unknown corrections to all the exterior 

orientation parameters (i.e. Aw^. , , Ak^. , AX^, aY^ and AZ^); and 

5. is a matrix of the unknown corrections to the coordinates of 
J 

point j, i.e. AX., AY. and aZ-, The following expressions for the 

J J J 

variance-covariance matrices were originally developed by Brown and 
have been used in common practice (2), (4); 

N"\ (2.3) 

and 

- 4 Nj'' * 

This method of the propagation of variance was a direct extension 
of the technique which was developed for direct linear least-squares 
regression, which has been treated extensively in textbooks on statistics 
(1), (7). Although many authors (8, 9) have discussed the problems 
of iterative least-squares adjustment in geodetic and photogrammetric 
literature, there has been no critical examination into the validity 
of Brown's formulation. 

This investigation showed that Brown's formulation is theoretically 
valid in direct least-squares solution, in which the original observation 
equations are linear functions of the unknowns. In iterative least- 
squares solution, as required by most photogrammetric problems, this 
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formulation is valid only if the unknown corrections (6 and 6^ in 
equation (2.1) and (2,2) respectively) converge to zero. In Section 
2.1, the theoretical basis of this formulation for a direct least- 
squares solution is examined; and in Section 2.2, the condition under 
which the formulation becomes invalid for an iterative solution will 
be discussed. 

2.1 Direct Linear Least-Squares Solutions 

2.1.1 A Linear Regression Model 

Consider a simple problem in linear regression. Figure 1 pre- 
sents the result of an experiment in which the value of a variable y 
was measured at various values of x. Let it be assumed that, in the 
absence of any error in the measurements, the variable y should be a 
linear function of x; i.e. 

y = a^ + a^x (2.5) 

Furthermore, let it be assumed that the measurement of the variable x 
can be made exactly. Thus, the only explanation for not having all 
the points in Fig. 1 lying in a straight line is the presence of 
random errors in the measurement of the variable y. Thus, for point 
i in Fig. 1, the following error equation can be generated; 

j'l «-jA^ 

where v^. is the random error in the measured value y^, Rearranging 
terms, the above equation may be expressed as follows: 

V, . + a,x,. - y, (2.6) 


[y,] 


In matrix notation. 


[V,] t [1 XjJ 


I.e, 


Vi ^ B) A = Y. 


(2.7) 
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Fig. 1. A Linear Regression Problem 
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With one equation generated by each measurement of y, m points will 
generate m equations which may be combined together as follows: 


'v,' 




''i 

>'2 

+ 

“2 

A = 

''2 

* 
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* 



m 


- 


fl 



“n, 


E 

1 


i. e, V + B-A = Y (2.8) 

It can be easily derived that the corresponding normal equations will 
be 

(B%) • A = B% (2.9) 

where W is the weight matrix for the measured values of y and is de- 
fined as follows: 

W = (2,10.) 

where is the standard error of unit weight; oy is the covariance- 
covariance matrix of the measured values of y. 

Let N = B^WB, and C = B% (2.11) 

The most probable value of the coefficients a^ and a^ can then be com- 
puted as 

A = N‘^C. (2.12) 

It can be theoretically proved that, for this solution under the given 
set of assumptions, the variance-covariance matrix for the computed 
coefficients may be simply expressed as 
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(2.13) 


2.1.2 Statistical Definition of Variance and Covariance 

Before attempting to derive Eq. (2.13), the statistical definitions 
of mathematical expectation, variance and covariance are first re- 
viewed (7). Let X be a random variable which can take on the values 
of Xp Xg, x^, ... Xj^ with the probability p^ P 2 > P3» ••• Pp 
respectively. Then, the expected value of X, denoted as E(X), is de- 
fined as 


n 

E(X) = 2 x.-p^ (2.14) 

i=l ^ ^ 


The expected value of X is also called the population mean of X. If 
X is a continuous random variable which has a probability distribution 
function f(x), then the expected value is defined as 


E(X) = 


CO 

xf(x)dx 

4 

00 


(2.15) 


It can be easily proved that if b is a constant scaler, then E(bX) = bE(X). 
The variance of the variable X is defined as follows: 

a^= E {(X-E(X))^} . (2.16) 

The covariance between two random variables X and Y is defined 
as follows: 

oj^y ={E (X-E(X)) (y-E(Y))} (2.17) 

Let Z be a column matrix of random variables Z-j , Zg, ... Z^^^, i.e. 
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Then, it can be easily derived from the above fundamental definition 
that the variance-covariance matrix for Z may be expressed as follows: 

= E{(Z-E(Z)) (Z-EU))"^} (2.18) 

where is an (mxm) matrix. 

2.1.3 Derivation of the Variance-Covariance Matrix 

The variance-covariance matrix of the unknown matrix A in Eq. 

(2.12) can thus be expressed as follows: 

= E{(A-E(A)) (A-E(A))^} (2.19) 

Substituting Eq. (2.11) and (2.10) into the above expression yields 

= E{(N~^C - E(N"^C))CN~^C-E(N”^C))'*’} C2.20) 

But E(N"^C) = E(N'^B^WY), Since Y is the only matrix consisting of 
variable parameters, 

E(N“^C) = ECN'^B^) = U'^B%(Y) 

Therefore, Eq. (2.20) may be expressed as 

= E{(N"^B% - N"^B%(Y))(N"^B%-N'’^B%(Y))^} 

= E{ (N'^ b'^W(Y-E(Y) ) (Y-E(Y))’’WBN'^ } 

= N'^ B%{ (Y-E(Y) ) (Y-E(Y) )^}WBN"^ 

By definition, 

(jy = E{(Y-E(Y))(Y-E(Y))'^} 

From Eq. (2.10), = a^W'\ 


Therefore, 
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= OqN'^NN"^ 

2 -1 

Hence a. = a„N ; 

A 0 

which is Eq. (2.13), 

In summary, it can be stated that in a given least-squares ad- 
justment problem, if the observation model may be expressed as 

V + B-A = Y; 

I 

and if the only random variables included in the matrix Y are the 
measured parameters; then the variance-covariance matrix of the un- 
known matrix may be computed as 



It will be shown in the Section 2.3 that the second condition is not 
necessarily satisfied in an iterative least^^squares solution, 

2.2 Iterative Least-Squares Solutions 

2.2.1 Simultaneous Solution for Aerotriangulation 

When the observation equations are non-linear functions, an 
iterative least-squares solution is generally required. As an example, 
consider the solution model for the simultaneous adjustment of photo- 
grammetric blocks. The mathematical formulation for this solution 
method has been well documented in the literature and will therefore 
only be presented here in summarized form. 

The solution is based on the following pair of col linearity 
equations which relate the image coordinates (x.., y. .) to the ground 

I J I J 

coordinates (X., Y- and Z.); 

J J J 
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aind 




F = V- ‘-y + 


[mi3(Xj-X^) + f"23(Yj.-Y^) + 

f[m^2(Xj-X^) + m22(Yj-Y^) + 
[iti^3(Xj-X^) + 


0 


( 2 . 21 ) 


0 


where Xp and yp are the image coordinates of the principal point j 
f is the focal length; X^, Y^ and are ground coordinates of the ith 
camera position during the exposure; and the niji's are functions of 

* J 

the three rotation parameters (uj^, <j>^, k^) of the camera. 

By first«^order Newton approximation, the collinearity equations 
may be linearized as follows: 



( 2 . 22 ) 


where 



n 


'ij 


■ ■'• ij P - - 


[»,3(X/-Xf ) . «-23(Y“-Y-) * .33(Z^-Zr)] 


0 vCO^ 


7O -7CO1 


and 


ij 


* [” 13 ^?-*“) " "' 23 ( 1 -''?°) * ” 33 ( 1 - 1 °)) 


(2.23) 


where X°, Y^> z9 are approximate ground coordinates of point j, 

J V w 

X^°, y 9°, Z^° are approximate ground coordinates of the ith 
camera position 


( — etc. are partial derivatives of Eq. (2.21) com^ 

3 yC "i 

‘'i puted using the approximate values of all the 

parameters. 

Equation 2.22 may be expressed in matrix notation as 



+ B. .A B. -A- 
U IJ J 



(2.24) 


where A is a matrix of the corrections to exterior orientation parameters 
(ax9, ay9, azJ, Aoj. , A^., AkJ of all camera stations; and *A*. is the 

1 1 1 I 1 1 J 

matrix of the corrections to the ground coordinates of point j. 

The complete set of observation equations for all measured image 
coordinates will take the following matrix form: 

« • « • 

V BA + Ba - e; 


which may also be simply written as 
V Ba = e 


(2.25) 


An observation equation may also be written for each coordinate of a 
known control point and for each known exterior orientation elements. 
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For example, the observation equations for a known control point with 
measured coordinates (X?°, z9°) are as follows: 

\i J 


V 



V 


AX. 


x° 


^^3 


y9 


AZ. 

J 





(2.26) 


In matrix notation: 




(2.27) 


In general, the complete set of control data may be represented by the 
following equation: 


V - 



(2.28) 


Therefore, combining Eq. (2,25) and (2.26), the complete mathematical 
model is as follows: 

V + BA * e 

•• •• 

V - A = e 


which may be simply expressed as 
V + Ba = e. 


The normal equations will be 
(B%)a = B^We 


(2.29) 


(2.30) 


where W is the weight matrix for the observations. 

The least-squares solution for A is then obtained by the solution 
of the above normal equations; i.e. 


A = 


(2.31) 
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where N = (B^WB) and C = B^We. 

However, because of the first-order approximation used in the 
linearization of the collinearity equations, the least-square solution 
must be iterated until the corrections (a) become negligibly small. 

It will be shown in the next section that the convergence of the 
A-matrix is the necessary condition for the expression 
0^ = (B%)"^ to be valid. 

2.2.2 Derivation of the Variance-Covariance Matrix 

The normal equation in Eq. (2,30) is in exactly the same form 
as the normal equation in Eq. (2.9), Therefore, the derivation pro- 
cedure should be identical to that presented in Section 2.1.3 until 
the following expression is reached: 

0^ = N’^B^WE{(e-E(£))(8-E(e))^}WBN"^ 

Since, by definition, E{ (e-E(e))(e-E(e))^} = , 

a = N'^B^Wa WBN'\ (2.32) 

A £ 

Two conditions of convergence of the least-squares solution are 
next considered separately, 
a) Convergence with a approaches zero. 

As A becomes zero, Eq. (2.29)becomes 

V = e. 

That is, the residual terms e and e in Eq, (2.23) and the residual 

ij -^ij 

terms e. in Eq. (2.27) are caused only by errors in the measured 
parameters. Therefore, the approximation parameters in these observation 
equations have become constant parameters and the only random variables 
in the residual terms are the measured parameters. Let denotes 
the variance-covariance matrix of all the observations. Then it is 
obvious that under this condition of convergence, the following relation- 
ship is true: 
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a = (2.33) 

e 

Since by definition, W = 

a = (7°° = al W'' (2.341 

e 0 

By substituting Eq. (2.34) into Eq. (2.32), it can be easily shown 
that again 



Let denote the matrix of all the approximate parameters, and X 
denotes the matrix of all the unknown parameters; i.e. 

X = X° + A 

Then, since X° has stabilized to become a constant term, the variance- 
covariance matrix of the computed values of the unknowns is equal to 
the variance-covariance matrix of the corrections; i.e. 

X A 

b) Solution stabilized but a fails to approach zero. 

In a weakly conditioned photogrammetric solution, such as that 
caused by poor geometry or low-accuracy controls, some or all of the 
correction parameters may never approach zero. After a certain number 
of iterations, these correction parameters may simply oscillate between 
certain boundary limits from iteration to iteration. If the iteration 
procedure is allowed to continue, such a solution will eventually begin 
to diverge rapidly. 

In such a solution, the residual terms in the observation equations 
will include measurement errors, approximation errors, and errors 
introduced by the linearization of the observation equations. It 
is no longer true that the variance-covariance matrix of these residual 
terms is equal to the variance-covariance matrix of the measurements; 

1 .e. 
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, 00 
a f a 
e 

2 -1 

Consequently, ^ N . 

There is no statistical method available for directly computing 
from the solution. 
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3. ACCURACY ANALYSIS IN LUNAR PHOTOTRIANGULATION 
3.1 Computer Program SAPGO-A 

A computer program was developed for performing complete error 
analysis on lunar phototriangulation. It is based on the principle of 
propagation of variance, as described In Chapter 2. The program may 
be used for either the triangulation of real photogrammetric data, or 
for simulation study using fictitious data. It can accept as input 
the following types of data: 

1) photographic image coordinates with a specified variance- 
covariance matrix for all points; 

2) ground controls with known rectangular ground coordinates 

(X., Y. and Z-) and the corresponding variance-covariance matrix 

\l J J 

for each control point; 

3) distances, horizontal angles and azimuths measured on the lunar 
surface with the corresponding standard error for each measure- 
ment; 

4) the position of any exposure station in rectangular coordinates 
(X^, Y^ and Z^) with its corresponding variance-covariance matrix; 

5) the orientation of the camera axis defined by the three rotation 
angles to, (p, and k at each camera station, together with the 
corresponding variance-covariance matrix; 

6) the straight-line distance between a camera station and any 
ground point, together with its standard error; and 

7) the known ground coordinates of a set of check points. 

The output of the program includes the following: 

1) the correction parameters at the end of each iteration; 

2) the computed standard error of unit weight at the end of each 
iteration; 

3) the final adjusted values of the exterior orientation parameters 

and their estimated standard errors (X^ + ay*', Y^ + Z*: + 

1 - X. 1 - y^- 1 " Z. 

“i ± <^i and ic + a ); 

1 ^1 i 
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4) the final adjusted ground coordinates of all the pass points 

and their standard errors; and 

5) the true errors at the check points. 

This program was a modified version of the SAPGO-A program which 
was originally designed for aerotriangulation in earth application. 

SAPGO is the acronym for simultaneous s<lJustment of £hotogrammetric and 
geodetic observations. The mathematical formulation of the solution 
has been published in (5). The modified version of SAPGO-A is equally 
applicable for both lunar and earth application, and a full documenta- 
tion of the program is being prepared (6). 

3.2 A Simulation Model 

A simulation was conducted to achieve two objectives: 1) to test 
the reliability of the method of accuracy analysis in the computer 
program SAPGO-A; and 2) to evaluate the contribution of the various 
types of control data on the accuracy of lunar phototriangulation. 

Fictitious data was generated by a computer program called 
DTAGEN, which was developed at the University of Illinois. For a 
given flight configuration specified by the camera focal length, flight 
altitude, percentage overlap, origin of the reference coordinate system, 
and position of the first exposure; program DTAGEN generates all the 
necessary photo image coordinates, exterior orientation parameters, 
and ground coordinates of all corresponding image points. To simulate 
the measurement errors in practice, these parameters were perturbed 
according to some specified standard errors. Both the perturbed and 
the true values of all of these parameters were provided by the program. 

The simulation model in this study consisted of a strip of 11 
photos having the same flight configuration as the 76-mm. metric 
photography from the Apollo 15 mission. The flight attitude was 
110 Km and the longitudinal overlap was 60%. There were 25 pass points 
arranged in a 5x5 rectangular pattern in each photo, and the photo 
coordinates were perturbed with a standard error of +^5 ym. 

Three types of control data were generated: 1) laser altimeter 
measurement; 2) attitude data , <t>^ , k^.) for the camera; and 3) 
tracking data. The laser altimeter measurements were the straight-line 
distances between the camera stations and the pass points located at 
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that center of the photos. These distances were perturbed with a 
standard error of +3 m. The camera attitude data were perturbed with 
the following standard errors; a ~ +10", a. = +20", and a = +10", 

OJ — 9 — K ~ 

These accuracy levels have been found to be realistic from the analysis 
of actual data from the Apollo 15 mission (7). 

In reality, the tracking data for each orbit consisted of a 
position vector in the epoch, velocity vector of the orbit, and time 
of exposure with respect to the epoch. These data can be used to 
constraint the position of the exposure stations. Since the objective 
of this study was to determine the internal accuracy of the photo- 
grammetric solution rather than the absolute positioning accuracy, 
the exterior orientation parameters (X , Y , Z , w, ((>, and k) of 
the first photo in the strip was held fixed and assumed to be error 
free. 

Helmering (7) reported that, from the analysis of three strips 
of Apollo 15 photography, the attainable accuracy for the velocity vector 
was about +0.1 meter per second while the timing accuracy was +1 
millisecond. In this study, simulation runs were performed using three 
different levels of accuracy in the velocity vector: +0.1, +0.5 and 
£1.0 meter per second. The relative timing accuracy was assumed to be 
+1 millisecond. At these accuracy levels, the corresponding uncertainty 
in the coordinates of the exposure stations corresponds to those listed 
in Table 1. All three coordinates (X^, Y^, of each exposure 
station were assumed to have the same standard error. 

3.3 Experimental Verification on the Propagation of the Variance 
3,3.1 Well -control led photography 

Figures 2, 3, 4 and 5 show the results of a series of simulation 
using various combination of control data, A tracking accuracy of 
+ 0.5 m/sec. for the velocity vector was assumed for all cases; i.e. 
the exposure station coordinates were weighted according to the standard 
errors listed in column 2 of Table 1. The following four combinations 
of controls were used: 



Table 1. 


Standard Errors of 
for Three Levels of 




1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


0 

+4 m. 
6 

7 

8 
9 

10 

11 

12 

12 

13 
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osure Station Coordinates 


= +0.5 m./sec. 


~ m./sec. 


0 

0 

+20 m. 

+40 

29 

57 

35 

70 

40 

81 

45 

90 

50 

99 

54 

107 

57 

114 

61 

121 

64 

128 
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1) tracking data alone {Fig. 2); 

2) tracking data and altimeter measurements (Fig. 3); 

3) tracking data and attitude measurements (Fig. 4); and 

4) tracking, altimeter and attitude data (Fig. 5). 

Five independent strips of 11 photos were triangulated. In these figures, 
the true errors in the adjusted coordinates of the pass points located 
near the center of each photo are plotted against the estimated RMS 
error computed from the propagation of variance. With the exception 
of a few data points, all the true errors felt within the three sigma 
value of the computed RMS errors. Thus, the RMS errors present a 
realistic evaluation of the accuracy of the computed coordinates. 

Tables 2, 3 and 4 list the estimated RMS errors versus the true 
errors in the computed coordinates of the exposure stations for five 
cases in which tracking, altimeter and attitude data were all used as 
controls. Also listed in these tables are the magnitude of the correction 
parameters in the last iteration of the iterative least-squares solution. 
Except in case 5 of Table 4, all the corrections were less than the estimated 
RMS errors of the parameters. 

3.3.2 Poorly Controlled Photography - Cantilever Extension 

Figures 6, 7 and 8 show the true errors in the computed coordinates 
of the center pass points for a cantilever strip. The exterior orienta- 
tion parameters of both the first and second photo in the strip were held 
fixed in the solution. The strip was not controlled by any other control 
data. Five independent strips were triangulated, and the true errors were 
plotted in these figures against the three-sigma value of the estimated 
RMS errors. 

These figures show rapid accumulation of systematic errors in the 
phototriangulation process. Furthermore, the systematic pattern varies 
from strip to strip. The computed RMS errors failed to detect these 
systematic errors and clearly underestimated the inaccuracy of the solution. 

Tables 5, 6 and 7 list the true errors in the computed coordinates 
of the exposure stations for the same cases. The magnitude of the 
corrections in the last iteration are also listed. It was observed 
that beyond station No. 7, the corrections generally exceeded the 
estimated RMS error of these parameters. 
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Table 2. Expected RMS Error vs. True Error in X-Coordinate of 
Exposure Stations (Tracking + Altimeter + Attitude 
Controls With = +0.5 m./sec.) 


Station 

Expected RMS 

True 

Error 

in Case 


Last 1 

Correction 

in 

Case 


Errors 
{+ meters) 

1 

2 3 4 

(meters) 

5 

1 

2 

3 

4 

5 

1 

- 

- 

- 

- 

- 


- 

- 

- 

- 


2 

10 

-5 

-10 

4 

26 

-6 

-3 

-3 

5 

-6 

1 

3 

11 

+3 

-13 

1 

6 

6 

-6 

-3 

6 

-6 

2 

4 

12 

25 

20 

18 

3 

-10 

-9 

-2 

5 

-4 

3 

5 

13 

-1 

-10 

-19 

-4 

-6 

-10 

2 

7 

-8 

2 

6 

13 

-2 

-32 

-28 

2 

14 

-10 

6 

7 

-2 

7 

7 

14 

-7 

-17 

-31 

2 

21 

-5 

-2 

4 

-7 

5 

8 

15 

3 

-14 

-39 

7 

-4 

-9 

-6 

8 

-6 

5 

9 

15 

-16 

-25 

-36 

9 

-25 

-n 

-7 

8 

-9 

0 

10 

16 

-16 

-11 

-38 

13 

-33 

-12 

-3 

10 

-9 

1 

11 ■ 

18 

13 

0 

-36 

17 

-8 

-9 

-2 

5 

-7 

0 
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Table 3. Expected RMS Error vs True Error in Y-Coordinates of Exposure 
Stations (Tracking + Altimeter + Attitude Controls With 

~ iO.5 m./sec, )' 

Expected Error in Case Last Correction in Case 

Station RMS 



Error 
{+ meters) 

1 

2 

3 

(meters ) 

4 

5 

1 

2 

3 

4 

5 

1 

- 

- 

- 

- 

- 

- 

- 

- 




2 

7 

-14 

4 

-16 

- 5 

+ 3 

-2 

3 

0 

1.0 

1 

3 

7 

11 

13 

- 9 

- 3 

- 5 

-4 

2 

0 

0 

1 

4 

8 

4 

11 

4 

-15 

-13 

-3 

4 

0 

0 

1 

5 

8 

-17 

5 

- 9 

-24 

- 4 

-5 

3 

0 

-1 

0 

6 

8 

8 

4 

- 4 

- 1 

0 

-5 

5 

0 

-1 

2 

7 

9 

6 

19 

- 6 

- 8 

- 8 

-4 

5 

0 

-1 

3 

8 

9 

12 

27 

- 5 

- 9 

-16 

-5 

8 

0 

-1 

3 

9 

10 

1 

19 

- 1 

-17 

11 

-6 

8 

0 

0 

3 

10 

n 

2 

16 

16 

-12 

- 8 

-5 

10 

1 

-1 

5 

11 

11 

16 

17 

10 

-15 

4 

-5 

9 

0 

-2 

4 
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Table 4. Expected RMS Error vs True Error in Z-Coord1nates of Exposure 
Stations (fracking + Altimeter + Attitude Controls With 
ay = +0.5 m./sec. ) 


Expected 
Station RMS 
Error 
(+ meters ) 


True Error in Case 
12 3 4 


Last Correction in Case 
5 1 2 3 4 5 


1 


2 

6 

10 

3 

8 

8 

4 

9 

13 

5 

9 

23 

6 

10 

23 

7 

n 

12 

8 

n 

15 

9 

12 

11 

10 

14 

11 

n 

16 

16 


(meters ) 


- 


- 

W 

- 8 

2 

- 4 

9 

-20 

- 6 

2 

0 

-13 

- 4 

0 

8 

-14 

- 1 

0 

11 

-15 

- 7 

- 5 

8 

-19 

9 

- 6 

13 

-19 

-18 

1 

5 

-12 

-14 

-14 

6 

- 6 

-18 

-14 

7 

- 8 

-21 

-23 

11 


- 2 

- 

0 

2 

1 

- 4 

- 1 

0 

3 

- 3 

- 4 

- 2 

-2 

3 

- 5 

- 3 

- 5 

-2 

5 

- 8 

- 3 

- 7 

^3 

5 

-13 

- 5 

-10 

-2 

4 

-17 

- 7 

-10 

-2 

4 

-20 

- 9 

- 8 

-4 

5 

-22 

- 8 

- 6 

-4 

8 

-25 

-12 

- 5 

-6 

10 

-24 
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Table 5. Expected RMS Error vs. True Error in X-Coordinates of Exposure 
Stations (Cantilever Extension) 

True Error in Case Last Correction in Case 


Station 

Expected 

RMS 

Error 

1 

2 3 

( meters ) 

4 

5 

1 

2 

3 

4 

5 

1 

2 

3 

31 

109 

70 

-19 

- 18 

- 5 

16 

-19 

-15 

47 

- 1 

4 

37 

100 

79 

61 

-137 

2 

16 

-34 

26 

17 

- 6 

5 

35 

74 

-188 

116 

-184 

12 

7 

-35 

1 

-19 

-11 

6 

35 

24 

-356 

171 

-190 

33 

16 

-13 

-32 

9 

- 5 

7 

36 

-105 

-467 

248 

-231 

93 

33 

- 7 

-20 

-19 

- 2 

8 

37 

-168 

-553 

244 

-325 

137 

59 

33 

2 

-15 

0 

9 

37 

-299 

-573 

169 

-454 

192 

71 

79 

19 

-43 

- 3 

10 

43 

-670 

-637 

111 

-590 

216 

74 

102 

68 

-83 

- 6 

n 

59 

-360 

-667 

- 4 

-753 

219 

57 

142 

77 

-103 

-16 
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Table 6. Expected RMS Error vs. True Error in Y-Coordinates of Exposure 
Stations (Cantilever Extension) 


Expected True Error in Case 

Station RMS 

Error 12345 

(meters ) 


1 

« 

- 


- 

- 

- 

- 

- 

- 

- 

- 

2 

3 

20 

42 

-15 

-17 

- 57 

1 

4 

-16 

- 6 

21 

- 5 

4 

26 

- 8 

-34 

-28 

- 92 

27 

5 

.27 

15 

24 

-10 

5 

30 

- 92 

-26 

27 

-127 

61 

- 4 

.17 

6 

17 

-17 

6 

29 

-139 

-35 

10 

-131 

84 

-16 

-11 

-21 

24 

-26 

7 

30 

.187 

-46 

- 12 

-127 

107 

-20 

-12 

-24 

36 

-34 

8 

30 

.233 

-13 

- 40 

-no 

133 

-30 

- 1 

-34 

59 

-50 

9 

32 

.277 

29 

- 58 

-113 

160 

-44 

19 

-25 

38 

-60 

10 

37 

.391 

77 

- 52 

.158 

195 

-52 

4 

-37 

65 

-81 

11 

46 

-542 

122 

.103 

-153 

223 

-37 

58 

-70 

35 

-99 


Last Correction in Case 
1 2 3 4 5 
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Table 7. Expected RMS Error vs. True Error in Z-Coordinate of Exposure 
Stations (Cantilever Extension) 


Station 


Expected 

RMS 

Error 


True Error in Case Last Correction in Case 

2 3 4 5 1 2 3 4 


5 


(meters ) 


1 

- 

- 

2 

- 

- 

3 

12 

- 36 

4 

19 

- 71 

5 

22 

- 85 

6 

23 

-104 

7 

23 

-114 

8 

23 

-101 

9 

26 

77 

10 

37 

- 31 

11 

55 

- 7 


21 

27 

-27 

11 

29 

1 

23 

17 

15 

87 

- 1 

7 

178 

93 

- 2 

285 

- 91 

-34 

381 

-141 

-17 

499 

-168 

- 7 

613 

-149 

27 


11 

- 5 

12 

34 

- 5 

21 

58 

0 

27 

85 

2 

34 

114 

- 4 

38 

125 

-14 

28 

124 

-26 

9 

128 

-36 

- 9 

140 

-39 

-55 


4 

-15 

5 

-17 

-16 

11 

-17 

- 2 

20 

- 3 

- 9 

31 

8 

- 2 

36 

22 

4 

39 

36 

18 

40 

49 

42 

40 

65 

93 

45 
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3.3.3 Conclusions from experiments 

Based on these experimental results, the following conclusions 
can be drawn: 

1. In phototri angulation problems where only low-accuracy controls 
are available, the correction parameters generally will not 
reach zero in the iterative least-squares solution. However, 
the correction parameters should converge to a value less than 
the estimated RMS errors computed from the propagation of 
variance. 

2. Generally, if controls are available throughout the strip or 
block (even though they are of low accuracy) and if the correction 
parameters converge to a value less than the computed RMS values, 
then the computed RMS values should be a good estimator for the 
accuracy of the computed parameters. 

3. The RMS errors computed by the propagation of variance cannot 
detect the rapid accumulation of systematic effects from random 
errors. It is well known that the double-summation effect of 
random errors produce systematic errors in phototriangulation. 

4. Generally, when the corrections in the last iteration exceed 
the computed RMS errors, then the RMS errors are not reliable 
as estimator of adjustment accuracy. 

It is recommended that in simulation studies, true errors of 
the adjusted parameters should always be computed to check on the 
computed RMS error For problems in which the controls- are sparsely 
distributed, several independent simulation cases should be performed 
to provide a check between the estimated RMS errors and the true errors. 

3.4 Effectiveness of Various Control Data 

The effectiveness of the three types of control data in improving 
the accuracy of lunar phototriangulation was studied by the propagation 
of variance. Three levels of tracking accuracy were assumed: jjD.l, 
jjD.5 and fl .0 m./sec. Altimeter measurements were assumed accurate 
to +^3 m. , and the attitude data were assumed to be accurate to +10 sec., 
+20 sec. and +10 sec. for u, (ji and k respectively. For each level of 
tracking accuracy, the following combinations of control data were con- 
sidered; 
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1) tracking alone; 

2) tracking and altimeter; 

3) tracking and attitude; and 

4) tracking, altimeter and attitude. 

Figures 9, 10 and 11 show the estimated RMS errors at the center 
pass points for the different combinations of controls at the three 
levels of tracking accuracy. 

The results showed that the altimeter measurements mainly improved 
the accuracy in the X-coordinates of the pass points; while the attitude 
measurements improved the accuracy of the Y-coordi nates. This result 
was anticipated because the altimeter measurement served primarily as 
scale control, while the attitude data served as direction control. 

Both types of control data were needed to achieve maximum accuracy in 
phototriangulation. Both the altimeter and attitude controls improve 
the elevation accuracy by about the same amount. 

The cases in which all three types of controls were used are of 
particular interest. It can be concluded from Figures 9, 10 and 11 
that for a single stereo model, the coordinates of the pass points 
can be determined with an RMS error of +7 m. in both X and Y and +15 m. 
in Z. At a tracking accuracy of +0.1 m./sec. for the velocity vector, 
there is practically no degradation of accuracy along the strip due 
to error propagation. At a tracking accuracy of +0,5 and £1,0 m./sec., 
the degradation in elevation accuracy along the strip is also negligible; 
but in both instances, the accuracy in the X- and Y- coordinates degenerates 
at the rate of about 1 m. and 0.5 m. per model respectively. 




Center Pass Point at Photo No. 
g. 9. Computed RMS Errors at Center Pass Points 
(Velocity Vector +^0,1 m./sec.) 





Computed RMS Errors in Coordinates (meters) 





Computed RMS Errors (meters) 
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Center Pass Point at Photo No. 



Fig. n. Computed RMS Errors at Center Pass Points 
(Velocity Vector +^1.0 m./sec.) 
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4. ACCURACY ANALYSIS IN ABSOLUTE ORIENTATION 
4.1 Compute Program THREED 

The accuracy of absolute orientation depends on the following 
factors: 1) accuracy of the model coordinates; 2) accuracy of the 
ground controls, 3) density and distribution of control points, 4) size 
of the area, and 5) scale of the stereo model. 

A computer program with complete error analysis capability has been 
developed for performing absolute orientation. The program is code 
named THREED. It may be used to perform any one of the following 
functions: 

1. To perform absolute orientation. 

The program takes as input the model coordinates of a set of model 
points and the ground coordinates of a group of control points. Both 
the model and the ground coordinates of the control points can be 
weighted individually according to their variance-covariance matrices. 
The program computes the seven transformation parameters (Xj, Yj, Zy, 

0 ), <}>, K and scale) and their estimated standard errors. The program 
also transforms the model coordinates of any pass point in the ground 
reference system and determines the standard errors of the transformed 
coordinates. 

2. To determine accuracy of absolute orientation by simulation. 

In simulation application, the program takes as input 1) the ground 
coordinates of a set of control points; 2) the variances of the ground 
control coordinates; 3) the variances of the model coordinates; and 
4) the scale of the model. The program then generates a set of model 
coordinates for the given ground points and perturbs them according 
to the specified accuracy of the model points. It then performs a 
regular absolute orientation solution and outputs the estimated 
standard errors of the seven absolute orientation parameters. 

3. To determine the uncertainty in the orientation of a surface 
defined by a set of triangulated pass points. 

The direct output of any phototriangulation solution is the ground 
coordinates of a set of pass points and their standard errors. Program 
THREED can be used to determine the uncertainty in the orientation 
(Xy, Yy, Zy, ti), K and scale) of the surface defined by the set of 
pass points. 



40 


A complete documentation of this program is being published 
separately (8), 


4,2 Mathematical Formulation 

The program THREED is based on the equations for three-dimensional 
conformal transformation which are as follows: 


- 





"j 


"'ll Hz "'13 




= X 

^21 ^22 ^23 



"j 


"'31 "’32 "’33 



. 






where 


Xj, Yj, 

X 

mil = 

mi 2 = 

"’is " 
^21 
®22 
®23 

^31 " 

Hz 

®33 


z. model coordinates of point j 

w 

Z. ground coordinates of point j 

\J 

Zj three translations 

scale factor 
cos <t> cos K 

cos u) sin K + sin w sin <)i cos k 

sin w sin k - cos w sin (j> cos k 

-cos <ti sin < 

cos ID cos K - sin 0 ) sin (p sin k 

sin (D cos K + cos w sin sin < 

sin (|> 

-sin to cos (|) 

cos ID cos (}) 


(4.1) 


After linearization by first-order approximation, Eq. (4.1) may be 
expressed as follows: 
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hi hz hs h4 hs he 

^21 ^22 *^23 *^24 ^25 ^26 

hi hz ha h 4 ^35 he 




dzz 

^^32 



,e. V. + B.A + B.A. 

J vj iJ J 


(4.2) 


The model coordinates of each control point generates one set of equations 
as in Eq. (4.2). For m control points, the complete set of observation 
equations may be expressed as follows: 



* *9 ^ 

i.e. V + BA + BA = e (4.3) 


In order to permit flexible weighting of the seven transformation 
parameters as well as the ground control coordinates, one set of 
observation equations is introduced for each. The set of observation 
equations for the transformation parameters are as follows: 




o 

X 

n 

1 

X 

—1 o 
o 

aYt 

= Y? 

voo 

■ h 

.Zt 


7O0 

Ao) 

. 0 

“ £jJ 

0 

0 

3 

1 
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V , - A<j) = ((>° - <J>°° 
q> 

- Ak = 

v" ^ AX = A° - X°° 

where X°, Y°, ... and x° are approximated parameters; and X°°, Y°°, 
and x°° are measured parameters. In matrix notation, this set of 
equations may be simply written as follows: 


V - A = e 
(7,1) (7,1) (7,1) 

The observation equations for the jth ground control point are as 
follows: 


X°-X°° 

J J 
yO yOO 


z°-z?° 


Again, the superscript (o) denotes approximation parameters and the 
superscript (oo) denotes measured parameters. Equation (4.5) may be 
simply written as 

V • ^ A . = e . 

J J J 

The complete set of observation equations for all control points are 
as follows: 


V - A = e 


(4.6) 
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Combining Eqs. (4.3), (4.4) and (4.6) yields the following observation 
model : 





• • 




P 

V 


B 

B 


A 


• 

e 

V 

+ 

-I 

0 



= 

e 

V 


0 

-I 


A 


* 

£ 



- 





- 


i.e. V + B4 = C 

The normal equations is then as follows; 

(4.7) 

(B%)a = B^WC 

(4.8) 


where W is the weight matrix for the observations. An iterative 
solution procedure must be followed. An initial set of approximate 
values is assigned to all unknown transformation parameters. The 
solution solved for the corrections and then apply the corrections to 
the approximations. The solution is iterated until a stable solution 
is reached. 

After the last iteration, the variance-covariance matrix (oj) 
of the computed transformation parameters is computed by the following 
expression; 

oj = al (B%)-^ (4,9) 

2 

where 0 ^ is the variance of unit weight. The ground coordinates of 
all other model points and the corresponding standard errors are 
computed by the following expressions: 

= X ^ "’zlJ'j ^ * Xt 

“ X * '" 22 ^} * * 'fj 

h ' X ^ "’23^J * "’33Xj> * h 
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o 2 = +" l 3 , Zj ) V \ 2 

J J 

2-2 2 

+ £-(sin^cos<)Xj + (sin^sinK)yj + cosifZj] x 

o o o 

+ E-coS(()Sintcx. - cosi}>cosicy.J x" a 
J J 

. J. .-2 2 . 2 .-2 2 . 2 -2 2 . 2 
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J J J * 
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J J 
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2,-2 2 
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J J J 
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J 
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J J 
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+ misX + ni23^ cJy_ + 11I33X + 07 

J J J ^ 

The m.. terms in the above equations are defined as in Eq. (4.1). 
^ J 


(4.10) 
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4.3 Accuracy Analysis by Simulation 

One major application of program THREED is the determination of 
absolute orientation accuracy by the method of simulation. Table 8 
lists the simulation results for twenty-five cases in which the accuracy 
of the ground controls ranged between 1 km. and 2 km. The standard 
errors of the model coordinates ranged between 2 and 200 meters. The 
dimension of the square area ranged between 2 and 2000 km., while 
the scale of the model varied from 0.4 to 1.2. In all cases, the area 
was assumed to have twenty-five control points arranged in a 5 by 5 
rectangular array. 

Based on the results of these twenty-five cases, an unsuccessful 
attempt was made to develop a prediction equation for each of the 
transformation parameters. The following general polynomials was 
used in a regression analysis: 

a = 3o + a^axY2 + + a3A + a^S + ago^y/ + ago^^^^ + a^A^ + agS^ 

(4.11) 

where 

OxYz standard errors of ground control coordinates 

a standard errors of model coordinates 

xyz 

A area 

S scale 


The coefficients a^, a-j , ... and ag were determined in a regression 
analysis using the data in Table 8. Although excellent fit was obtained 

, and 07 ), errors in the pre- 
ters (a , ff., and a ) generally 

CO (p K 

exceeded +^1°, 

Furthermore, Eq. (4.11) gives no consideration to the density and 
distribution of ground controls, which, are important factors governing 
the accuracy of absolute orientation. 

Program THREED should provide the best means for evaluating absolute 
orientation accuracy. 


for the translation parameters {a^ , ay 
dieted accuracy of the rotational Jarami 



Table g- Accuracy of Absolute Orientation 


V 7 

a 

Di mens i on 
n T 1 

Mode 1 




Standard Error 

of Absolute Orientation Parameters 

A , 1 , Z. 

(meters) 

^ » y » z 

(meters) 

^1 ^ u u a 1 ^ 

Area 

(km) 

Scale 

Scale 

(meters) 

(meters) 

(meters) 

(ii 

0 In 

“ 1 II 

1C 

0 1 

1000 

20 

200 

0.8 

.0013 


17A 

17A 

0- 8-55 

0- 6-57 

0- 7- 0 

1250 

6 

20 

0.6 

.012 

219 

219 

219 

1-A8-51 

I-2A-5A 

1-29-31 

1250 

6 

20 

1.0 

.020 

217 

217 

217 

1-A7-59 

1-23-53 

1-28-33 

1250 

6 

2000 

0.6 

.0001 

218 

218 

218 

0-01-07 

0- 0-52 

0- 0-52 

1250 

6 

2000 

1.0 

.0002 

216 

216 

216 

0-01-07 

0- 0-53 

0- 0-52 

1250 

63 

20 

0.6 

.012 

217 

217 

217 

1-A8-37 

1-2A-38 

1-29-07 

1250 

63 

20 

l.O 

.021 

218 

218 

218 

1-A8-52 

1-2A-56 

1-29-30 

1250 

63 

2000 

0.6 

.0001 

216 

216 

216 

0- 1-07 

0- 0-52 

0- 0-52 

1250 

63 

2000 

1.0 

.0002 

217 

217 

217 

0- 1-07 

0- 0-52 

0- 0-52 

1750 

6 

20 

0.6 

.017 

305 

305 

305 

2-31-01 

1-59-13 

2-05-OA 

1750 

6 

20 

l.O 

.027 

296 

296 

297 

2-25-31 

1-53-03 

1-58- 3 A 

1750 

6 

2000 

0.6 

.00017 

30 A 

30A 

30A 

0-01-33 

0-01-13 

0-01-13 

1750 

6 

2000 

1.0 1 

.00028 

295 

295 

295 

0-01-3A 

0-01-13 

0-01-13 

1750 

63 

20 

0.6 

.018 

305 

305 

305 

2-32-20 

1-59-56 

2-05-1 A 

1750 

63 

20 

1.0 

.029 

305 

305 

305 

2-31-28 

1-59-A7 

2-0A-A7 

1750 

63 

2000 

0.6 

.00017 

303 

303 

303 

0- •1-33 

0- 1-13 

0- 1-13 

1750 

63 

2000 

1.0 

.00028 

30A 

30 A 

305 

0- 1-3A 

0- 1-13 

0- 1-13 

1500 

20 

200 

0.^4 

.001 

261 

261 

261 

0-13-21 

0-10-26 

0-10-29 

1500 

20 

200 

1.2 

.003 

261 

261 

261 

0-13-22 

0-10-26 

0-10-29 

1500 

20 

2 

0,8 

.12 

275 

278 

27A 

10-21-52 

10-A8-23 

15-A8-A7 

1500 

20 

2000 

0.8 

,00002 

261 

261 

262 

0- 0- 8 

0- 0- 6 

0- 0- 6 

1500 

2 

200 

0.8 

.0017 

2AA 

256 

250 

0-1 1-5A 

O-O9-3A 

0-09-27 

1500 

200 

200 

0.8 

.0019 

260 

260 

260 

0-13-21 

0-10-26 

0-10-28 

1500 

20 

200 

0.8 

.0019 

261 

261 

261 

0-13-23 

0-10-27 

0-10-29 

2000 

20 

200 

0.8 

.0026 

3A8 

3A8 

3A9 

O-I7-A9 

0-13-56 

0-13-59 
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4.4 Analysis of Apollo 15 Data 

Program THREED was used to evaluate the mapping accuracy of the 
Apollo 15 photography. Figure 12 shows the overlapping coverage be- 
tween a stereoscopic model of the 76 mm. metric photography and three 
successive stereoscopic models of panoramic photography. 

The stereoscopic model of 76-rran. metric photography was absolutely 
oriented on the AS-ll-Bl using orientation parameters derived from 
tracking and stellar camera data. The ground coordinates of the 
model points in Fig. 12 were measured directly from the model in the 
AS-ll-Bl. Table 9 presents a listing of these ground coordinates. 

The three panoramic stereomodels were also individually oriented in 
the AS-ll-Bl. The model coordinates of the control points as well as 
about 20 other model points were recorded. To minimize the errors 
introduced by the geometry of the panoramic photography, only the 
central 17.5® sweep portion of each frame (+8.25® from center) was 
used. Tables 10, 11 and 12 list the model coordinates for the three 
panoramic stereo models. 

All the above AS-ll-Bl measurements were performed by Lockheed 
Electronics Company, at the NASA Johnson Space Center, Houston. 

Program THREED was used to perform absolute orientation of the 
panoramic models and to perform error analysis. After several trial 
solutions, the following sets of standard errors for both the model 
and the ground coordinates were found to give the best agreement 
between the computed standard error of unit weight and the predefined 
value ; 


Model Coordinates 

= + 0.010 mm. 
c?y = + 0.010 mm, 

= + 0.010 mm. 


Ground Coordinates 

= ;+ 1 5 m . 

Oy = i 10 m, 

0 Z = + 25 m. 


Table 13 lists the standard errors computed by program THREED for both 
the transformation parameters and the adjusted ground coordinates. 
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It can be concluded from this study that the ground coordinates 
derived from this procedure have the following relative positional 
accuracy within the 76-mm, stereo model: = 4J5 m., ay = +^10 m. 

and 0 ^ = +25 m. The orientation accuracy of the panoramic stereo 
models can be expected to be +3 minutes of arc in the ui and k ro- 
tations and +^7 min. in the i> rotation. 
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Fig, 12, Stereoscopic Coverage at One Metric Model and 
Three Panoramic Models 
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Table 9. Ground Coordinates Derived From the Metric Model 


tiFOU^r CUORDINATtS 


PUINT X Y 


b 

*36^97 , 0000 

-9l 17,0000 

232,0000 

9 

-29682,0000 

-5536,0000 

377,0000 

10 

-22551 . 0000 

-8501,0000 

422,0000 

21 

-30412,0000 

13031,0000 

372,0000 

22 

-22827,0000 

6572,0000 

203,0000 

23 

-35327,0000 

22215,0000 

1182.0000 

?A 

-22453,0000 

21101 ,0000 

295,0000 

11 

-16923,0000 

-4620,0000 

1006,0000 

12 

-9550, OOOO 

-9020,0000 

476.0000 

25 

-17108,0000 

20768,0000 

338,0000 

26 

-9752,0000 

18663.0000 

326,0000 

19 

-9076,0000 

6680,0000 

471 ,0000 

20 

-19326, 0000 

8644,0000 

133,0000 

18 

-2361,0000 

5845.0000 

305,0000 

27 

1806,0000 

15096,0000 

310,0000 

29 

10759,0000 

21624,0000 

386,0000 

15 

14080,0000 

-9575,0000 

904,0000 

M 

8007 , 0000 

-6116,0000 

360,0000 

1 7 

6701.0000 

5200,0000 

301,0000 

30 

17232,0000 

9647,0000 

347, 0000 
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Table 10. Model Coordinates from Panoramic Model 22/27 



M D D 

el COO 

R 0 I N A T 

L S 

1, CONTRUL 

POINTS 




numpep 

POINT 

X 

Y 

Z 

1 

b 

-3b.flbOU 

-99, 39ftU 

1.6360 

? 

9 

V.0b20 

-76.4540 

2,4760 

3 

10 

56, 3bOO 

-96,6570 

2,6570 


21 

6, 0690 

47,5630 

2,4370 

5 

?2 

5b. 3340 

3,6140 

1,4570 

6 

?3 

-2b. 5300 

109.1900 

7,6360 

7 

?4 

59.9610 

100.4220 

2.2830 


2, unknown 

POINTS 




NUMPEP 

POINT 

X 

Y 

7 

8 

126 

-36,0220 

-17,6710 

2,3150 

9 

116 

-2.4400 

-18.2550 

3,0560 

1 0 

110 

26,1510 

-17,7480 

2,3440 

1 1 

103 

6J. S470 

-17,3580 

1.4610 

12 

1C2 

62.5500 

19,5060 

1,2460 

13 

109 

32,3130 

21,0880 

2,4340 

14 

117 

-0, 9240 

20.4110 

2,6310 

15 

125 

-36.9070 

20,1300 

2.5700 

16 

124 

-39,3250 

49.6360 

2,2550 

17 

108 

29,2470 

51,8240 

2,3590 

18 

101 

61.7660 

53.9560 

2,2960 

19 

107 

32,2640 

90,2580 

2,7140 

20 

115 

-2,0260 

67.3500 

4,1060 

21 

123 

-35, 3630 

67,5960 

2.6990 

22 

122 

-36,4530 

1 15.8250 

7,9560 

23 

106 

26,5690 

1 18,6400 

2,7150 

24 

127 

-33, 7020 

-44.6800 

2,3020 

25 

119 

-3.5630 

-50.3710 

4,1170 

26 

111 

40.7050 

-46,0320 

6,0690 

27 

104 

65,3690 

-47,8490 

4,3590 

28 

105 

64, 3700 

-65,4150 

4,9000 

29 

112 

36,4960 

-63.3250 

2,6670 

30 

120 

-2,0950 

-69,6810 

2,7150 

31 

128 

-35,1640 

-77,9410 

2,0910 
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Table 11, Model Coordinates From Panoramic Model 



M 0 

n f L CO 

n P 0 I M 

1, CONTROL 

POINTS 



NUMPr'» 

POINT 

X 


1 

11 

-37.9920 

-74,6960 

2 

12 

10,6850 

-103,1930 

3 

25 

-37.4780 

96,1670 

4 

26 

11.5460 

82,7240 

5 

19 

15. 3690 

2,6830 

6 

20 

-53,?6?0 

16,7970 

7 

18 

59.7380 

-4,6610 


2, UNKNOWN 

POINTS 



NUMBER 

POINT 

X 


8 

130 

-43,3g00 

110,1490 

9 

136 

-8,4390 

112,9540 

10 

144 

25.6170 

116,501 0 

1 1 

152 

59,1410 

111,1180 

12 

154 

59,6470 

40,5940 

13 

146 

24,7540 

50,9030 

14 

138 

-6, 1570 

56,0330 

15 

131 

-43, 7580 

45,321 0 

16 

132 

-43,7830 

10,4920 

17 

139 

-6,0850 

17,9130 

18 

147 

31 , 5370 

21,2860 

19 

155 

58.6940 

10,7360 

20 

156 

59.9080 

-22,9870 

21 

148 

27,6720 

-30,3710 

22 

140 

-3,7810 

-21,7520 

23 

133 

-38,9030 

-22,6870 

24 

134 

-42,6650 

-49,6090 


20/25 


t s 


2 

•2. Hio 
•5»5940 
-6.7080 
-6.8140 
-5.6740 
-7.8810 
-6,9810 


2 

•6 ,6990 
-6,4160 
-6,2930 
-6,2430 
-6.5840 
-6,3250 
-6,6140 
-8,0040 
-7,8070 
-6,7400 
-6,3660 
-6,4890 
-6,6030 
-6,1090 
-6,6150 
-7,83?0 
-6,1500 
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Table 12, Model Coordinates From Panoramic Model 18/23 



M 0 

n F L CO 

0 P 0 I M A 

T L S 

1, CONTROL 

POINTS 




NUMBFP 

POINT 

X 

y 

Z 

1 

27 

-45, Jo60 

56,3960 

-0,7990 

? 

?v 

14,9640 

100,6650 

-0,0610 

3 

15 

3J.9300 

-107.6350 

3, 1330 

4 

14 

-6,4020 

-63,6530 

-0,2150 

5 

17 

-0,4790 

-8,4340 

-0.7040 

6 

30 

56.B220 

20,2310 

-0,3690 


2, UNKNOWN 

POINTS 




NUMBER 

PUINT 

X 

Y 

7 

7 

2 

49,5940 

3<) , 1 320 

3,4940 

8 

26 

-36.9190 

106,5860 

-0,4200 

9 

175 

-6. 0570 

-122,2780 

0.275C 

10 

182 

19,1600 

-117, 1890 

3,0960 

11 

190 

56,6780 

-115,5610 

4,2270 

12 

181 

19,5660 

-84,5990 

1,9580 

13 

174 

-10,9570 

-68.4420 

-0.3510 

14 

166 

-42,5250 

-66,6240 

-0,2290 

15 

172 

-10, .>560 

-21.7630 

-0,5630 

16 

173 

-13,2660 

-57,2540 

-0,3670 

17 

180 

20,1330 

-50,1990 

0,0560 

18 

186 

57. 1720 

-55,6450 

2.1830 

19 

169 

-9,6530 

72.5150 

-0,1700 

20 

177 

20,3240 

75.1330 

-0.1280 

21 

184 

54,6090 

76,3210 

3,1110 

22 

105 

58,5690 

47,8650 

2,9140 

23 

1 76 

26,0900 

52,0830 

0,0280 

24 

170 

-10,6650 

51.4980 

-0,2360 

25 

163 

-41,6560 

5,0400 

-0,6790 

26 

171 

-5, 2450 

16,6500 

-0,1260 

27 

186 

54,8420 

13,6500 

-0.2230 

20 

187 

53,0320 

-14.1090 

0,5070 

29 

199 

26,9790 

-23.7330 

-1.1620 

30 

166 

-11 ,2190 

100,6190 

-0,0410 

31 

183 

56,2960 

112.7730 

2,5650 

32 

217 

23,8970 

110,4950 

0,3390 
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Table 13. Accuracy in the Absolute 
Model s 

Parameters 

Predefined 
Computed 



a 

LO 


a 

K 

Adjusted Control 
ay Coordinates 

Transformed 

ay Model Coordinates 

No, of Control Points Used 
in Solution 


Orientation of the Panoramic 


Panoramic Model 


18/23 

20/25 

22/27 

jjO.Ol mm. 

jfO.Ol mm. 

+0.01 mm. 

0.012 mm. 

0.010 m. 

0.012 mm 

7 m. 

6 m. 

7 m. 

5 m. 

4 m. 

5 m. 

13 m. 

9 m. 

12 m. 

3' 53" 

3' 4" 

3' 2" 

8' 47" 

5' 36" 

6' 49" 

1' 54" 

r 25" 

r 32" 

10 m. 

9 m. 

10 m. 

8 m. 

6 m. 

8 m. 

25 m. 

18 m. 

23 m. 

13 m. 

9 m. 

11 m. 

10 m. 

7 m. 

8 m. 

32 m. 

22 m. 

26 m. 

6 

7 

7 
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5. SLOPE ACCURACY 

5,1 Error Propagation Formula 

Let D be the horizontal distance between points i and j which have 
elevation h^. and hj respectively. Then, the slope angle from point i 
to point j can be computed by the expression 


^ij = = tan"'' (^) (5.1) 

Hence, by the law of propagation of random errors, the standard error 
of the slope may be expressed as follows; 


Since 


ij 


2 = a 2 + „ 2 

^9Ah ^ Ah 


(5.2) 


3Ah 


Ah 2 -1 1 

(1 ) J 


and 


as. . 

TJ 

aD 


Cl ♦ (- 4)- 

u 


5. • 

IJ 


[1 + (#)V^ 


C2 '^Ah ] 


Ah^i 


(5.3) 


Assuming that ^h’ 

J 


then 


2-9 
“ih - 2“h- 


(5.4) 


Furthermore, since D = E(X^-XO + (Y.-Y-) ] 

j I j ^ 


2nl/2 


J J 1 1 


Assuming that 



56 


then 


1 .e. 


h 

1 _ n2r,n2i 2 

Cq - D {2D } 

2 _ 2 


= a 


X,Y» 


Substituting Eqs. (5.4) and (5.5) into Eq. (5.3) yields 


_ /o r-i . / h\2'i~l rl 2 . /Ah\2 2-il/2 


(5.6) 


in which a<- is given in radians. In general, Eq. (5,6) may be re- 
^i i 

written as follows: 


^S. . 
ij 


= c [1 t [1^ ^ (^j2 ^2jl/2 


^X,Y 


(5.7) 


where c = 1,414 for Oc in radians; 

^ij 

= 81.016 for Oc in degrees; 
ij 

= 4861 for Og in minutes; 
iJ 

5 

= 2.916x10 for in seconds. 

ij 


When D is much greater than Ah, the term for may be neglected resulting 

in the following simplified expression; 



D[1 + (^)^J 


( 6 . 8 ) 


5.2 Application on the Apollo 15 Data 

Table 14 lists the THREED output ground coordinates for all pass 
points in the panoramic model 22/27. Table 15 lists the RMS errors in 
slopes along two profiles in this model. In all instances, the term 

an^ ill Eq. (5,3) is negligibly small and Eq. (5.8) was used for 
D 

computing the RMS errors. The parameter was computed as the mean of 
the RMS elevation errors of the two end points of each line. 



Point 


126 

lie 

110 

103 
102 
109 
117 
125 
1P« 
loe 
101 
107 
115 
123 
122 
106 
127 

119 
111 

104 

105 
112 

120 

ipe 


Table 14. Qround Coordinates of Pass Points in Panoramic Model 22/27 


Y 

Z 

X 

RMS Errors 
Y 

(meters) 





-36*576.465 

3120. 151 

339.285 

-31540.133 

3146.194 

453.093 

-26952.520 

3297, 8o8 

348,501 

-21626,826 

3441,667 

221.624 

-21682.055 

6968,773 

184.357 

-26421,641 

9133. 398 

360,276 

-31390.719 

8952.133 

387.427 

-36602.652 

6823. 164 

376.983 

-37236,326 

13243.367 

326.652 

-26955.492 

13736,559 

347,189 

-22079,617 

14134.652 

340.022 

-26595.426 

19509.113 

398,657 

-31732.160 

18990.473 

605,118 

-36733.367 

16947,098 

391.567 

-36965. 1«0 

23179.227 

1178,684 

-27214,973 

23757, 7o7 

397, 051 

-36165.941 

-890,930 

339.056 

-31634.402 

-1672.051 

613.865 

-24996,781 

-1214.539 

912.851 

-21294.375 

-1127,785 

655.118 

-21356.891 

-6765.205 

738.173 

-25543,000 

- 6518.941 

431.012 

-31364,633 

-4565.129 

404.653 

-36305.215 

-5883.762 

309,035 


7.373 

5,600 

16.933 

7,151 

5.012 

12,687 

7.274 

5.341 

15,215 

7.626 

6,550 

22.828 

7,790 

6.508 

22.471 

7,343 

5, 47o 

15.918 

7,162 

5.028 

12.732 

7,413 

5.596 

16.855 

8,044 

6.132 

18,420 

7.946 

5.910 

16.557 

8.456 

6,967 

23.155 

9,363 

7. o76 

19,578 

9, 175 

6.666 

17.015 

9.352 

7 . o7 8 

20. I 8 O 

10.677 

6. i6l 

22.754 

10.642 

e. 04 1 

21.600 

7.797 

5,861 

16.937 

7.646 

5,575 

14.142 

8.195 

6. 184 

18.612 

8.523 

f .040 

24,095 

9.771 

7.898 

25.733 

9.194 

6,969 

20.037 

8.425 

6.081 

15.493 

8,675 

6,^2 

19,095 
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Point 

From 

No. 

To 

Table 15. 

Horizontal 
Distance (D) 

RMS Errors in Slope 

Elevation Slope 

Difference (Ah) 

* 

^h 

RMS Error 
in Slope 



(meters) 

(meters) 




(meters) 


122 

123 

4,238 

-787 

-10° 

31 ' 

+21 

+23' 

123 

117 

11,333 

- 5 

- 

0° 

2' 

16 

1 7' 

117 

no 

7,188 

- 39 

- 

0° 

19’ 

14 

+10' 

no 

111 

7,183 

+565 

+ 

4° 

30' 

17 

+11' 

111 

105 

5,637 

-175 

- 

r 

47* 

22 

+19' 

115 

117 

10,044 

-218 

- 

r 

15' 

15 

+ 7' 

117 

118 

5,806 

+ 66 

+ 

0° 

39' 

13 

+11 ■ 

118 

119 

4,821 

+161 

+ 

1° 

55' 

13 

+13' 

119 

120 

2,906 

-209 

- 

4° 

07' 

15 

+25' 


*a|^ was computed as the mean of the RMS elevation errors of the two 
end points. 
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6. CONCLUSIONS 

Theoretically, the method of propagation of variance is a reliable 
means of evaluating the accuracy of an iterative least-squares solution 
only if the correction parameters converge to zero. However, simulation 
results in lunar phototriangulation showed that the method is also 
applicable to iterative least-squares solutions for which the follow- 
ing two conditions are achieved: 1) the corrections in the last 
iteration are smaller than the estimated standard errors of these parameters; 
and 2) the photogrammetric solution is well controlled to prevent any 
rapid accumulation of systematic errors. 

Two computer programs have been successfully developed to use 
the method of propagation of variance to evaluate the relative and 
absolute accuracy of lunar phototriangulation. Programs SAPGO-A 
and THREEO form an accuracy analysis package which may be used to per- 
form the following functions in any phototriangulation project: 

1) to predict the internal accuracy of photogrammetric solution 
(SAPGO-A); 

2) to determine whether auxiliary data or ground controls of a 
specific accuracy can strengthen the solution (SAPGO-A); 

3) to predict the accuracy of the final solution, including standard 
errors of both the triangulated pass points and the photo orienta- 
tion parameters (SAPGO-A); 

4) to determine the uncertainty in the tilt of the terrain surface 
defined by the computed pass points (THREED); and 

5) to determine the absolute orientation accuracy of any stereo- 
model (THREED). 

A simple expression was also developed for computing the standard 
error of the slope as a function of the standard error of the elevation 
at the two end points. 

Simulation studies of the Apollo 15, 75-mm. photography showed 
that both the laser altimeter and the stellar camera measurements 
contributed significantly to improving the phototriangulation accuracy. 

A relative accuracy of +7 m. in the X and Y coordinates and +15 m. in 
the Z coordinates should be obtainable for the triangulated pass points 
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in a single stereo model even if the velocity vector has an RMS error 
of +1.0 m./sec. At a tracking accuracy of +0.1 m./sec. in the velocity 
vector, there was practically no degradation of accuracy due to error 
propagation in the phototriangulation process. At a tracking accuracy 
of +0,5 m. and +1.0 m. sec,, the degradation in elevation accuracy is 
also negligible; but the accuracy in the X and Y coordinates degenerated 
at the rate of about 1 m. and 0,5 m. per model respectively. 

Analysis of the AS-ll-Bl measurements provided by NASA showed that 
the relative positional accuracy of the control points derived from one 
model of 76-mm. metric photography were: = £15 m., oy = ±10 m. and 

©2 = ±25 m. These RMS errors are approximately twice as large as those 
predicted from the simulation studies. This finding should not be sur- 
prising, because the simulation studies assumed a rather small RMS 
error of +5 ym. for the image coordinates; whereas, the model coordinates 
actually measured from the AS-ll-Bl were found to have an RMS error of 
+10 ym. A lower image fidelity together with some unavoidable instrumental 
and human error in the AS-ll-Bl measurements would then easily account 
for the larger RMS errors in the computed ground coordinates. 

The control points derived from the 76-mm. photography could be 
used to orient the high-resolution models of panoramic photography with 
a relative accuracy of +3 minutes of arc in the a> and < rotations and 
±7 minutes of arc in the ^ rotation. Within the ±8.25® sweep area of 
the panoramic model, the relative position of any point in the model 
area could also be determined with an RMS error of ±15 m. , ±10 m. and 
+25 m. in the X, Y and Z coordinates respectively. 
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